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William Parton, Steven Running, and Brian Walker 

Background 

Simplified or toy models have been commonly used to test specific 
hypotheses about complex systems. During the 1990 Global Change 
Institute (GCI). we attempted to propose the structure of various toy 
models of different earth system components from which we can 
eventually test hypotheses relevant to global change. These toy mod- 
els have been pared down to incorporate only essential attributes of 
key components and processes. 

The value of these toy models is derived from their relatively trans- 
parent nature. That is, processes connecting various components are 
easily discerned, and responses of these models to perturbations to 
their parameter environment can be more readily determined. The 
development and testing of these simplified models will quicken the 
rate at which we are able to unravel specific interlinkages of the 
earth system. Thus, these toy models can be viewed as an Important 
step in building more complex and intricate earth models. 

Scientists presented several toy models in discussions during the 
1990 GCI. These were: ocean biogeochemistiy (W. Broecker), global 
nitrous oxide (N 2 0) flux (I. Fung), biosphere interaction with climate 
(F. Bretherton), ocean carbon (B. Moore), terrestrial biogeochemistiy 
(W. Parton), forest carbon and water flux (S. Running), and sea ice 
growth (A. Thorndike). 

Carbon Flow Model 

At the GCI, a conceptual framework was developed for a general- 
ized carbon flow model for the major terrestrial ecosystems of the 
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world. The model is a simplification of the Century model (Parton et 
al. f 1987, this volume) and the Forest-Biogeochemical (BGC ) model 
(Running and Coughlan, 1988). The C flow diagram (Figure 1) 
shows that live foliage, roots, and woody biomass are represented 
and that C flows from these compartments to the plant residue 
compartment. Plant residue then decomposes and forms soil 
organic matter. The soil organic matter is divided up into two com- 
partments (slow and passive), with the slow pool having an approx- 
imate turnover time of 20 to 50 years and the passive pool having a 
long turnover time (1000 to 3000 years). The actual turnover times 
of these pools depend on abiotic decomposition factors (soil temper- 
ature and water), with shorter times in warm, humid environments 
and longer in cold, dry locations. This structure is a simplification 
of the Century soil organic matter model where the slow compart- 
ment of soil organic matter is further divided up into two boxes 
(active and slow with approximate turnover times of 2 and 20 
years, respectively). The model also represents the dynamics of 
nutrients in the plants and soil. The nutrient flows are represented 
because either N or P limits plant production for most terrestrial 
ecosystems, and production can be constrained by the ratios of C.N 
or C:P in the various components. The nutrient flow model has the 
same structure as the C flow diagram and will represent the 
dynamics of the most limiting nutrient for plant production. Nitro- 
gen limits most terrestrial systems; however, P limits production for 
many tropical forest systems. 


Carbon Dioxide 




Figure 1. Conceptual framework for a global terrestrial C model. 
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Plant Production 

Plant production Is calculated using two different submodels. The 
model developed by Running and Coughlan (1988) is used to calculate 
net daily photosynthesis as a function of soil water content, solar radi- 
ation, air temperature, and leaf area index. The calculated net photo- 
synthesis rate is the maximum potential plant production rate for the 
day: the actual plant production rate will be lower than the maximum 
rate il there is insufficient available inorganic N or P. The C fixed each 
day is then allocated into the live shoot, root, and wood biomass The 
allocation of C is primarily a function of the type of ecosystem being 
represented; however, the allocation pattern is also a function of the 
annual rainfall, grazing level, and atmospheric CO z level. 

The second submodel uses a less mechanistic approach where 
the controls on potential plant production are monthly rainfall and 
soil temperature. This type of model is currently being used in the 
Century plant production model. The effect of rainfall on plant pro- 
duction (Figure 2a) is calculated as a function of the ratio of 
monthly rainfall to potential evapotranspiration (increases linearly 
as the ratio increases). The potential evapotranspiration rate is cal- 
culated as a function of the maximum and minimum air tempera- 
ture using the equations developed by Linacer (1977). Production 
1 , I l Cr ^ aSeS rapidl y as surface soil temperature increases from 5° to 
5 C ’ then increas es relatively slowly up to 20°C, followed by a 
ecline in production above this temperature (Figure 2b) For cer- 
tain specified plant communities adapted to cooler climates, a curve 
similar to the C3 curve of the Century model can be used, or for 
wanner climates the C4 curve can be specified. The potential plant 
pro uction rate is calculated by multiplying the temperature effect 
times the rainfall term times the maximum plant production rate for 
a month. The actual plant production will be less than the potential 
rate if sufficient nutrients are not available. The allocation of C into 
the d | ff erent plant parts is calculated using the same approach as in 
the high-resolution plant production model. 

The live wood, foliage, and root biomass grow as a result of inputs 
of carbon from the plant production model. The different plant parts 
ave different lignin contents and nutrient concentrations, which 
are a function of the plant community. We have tentatively identified 
eight ecosystems that need to be represented in the worldwide appli- 
cation of the model: * * 

• Desert 

• Arid and semiarid savannas and shrublands 

• Humid savannas and woodlands 
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mixture of C3/C4 plants. 


• Grasslands 

• Temperate deciduous forests 

• Temperate evergreen forests 

• Arctic and subarctic tundra 

• Tropical rainforests. 
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The live roots die at a constant rate that is a function of the ecosys- 
tem-specific annual fine root turnover rate. The large wood biomass 
includes the boles, large branches, and coarse roots. The model con- 
siders the effects of harvest and natural death of the live wood. Har- 
vest of wood is an external driving variable, and wood death and 
decomposition are simulated. Wood death rates are ecosystem-spe- 
cific and a function of the amount of wood biomass (higher for older 
and larger trees). Leaf death rates vary as function of the vegetation 
type, air temperature, and moisture stress. The ecosystem-specific 
patterns for leaf death (e.g., deciduous vs. coniferous) and the con- 
trols of air temperature and moisture stress Eire represented. Burning 
redistributes carbon and nutrients from live foliage and wood to litter 
and the atmosphere (as oxides of nitrogen and C0 2 ). The redistribu- 
tion pattern needs to be modeled externally by a separate fire model 
with the results being fed back into the carbon model. 


Decomposition and Nutrient Cycling 

Dead leaf and fine root material flow into the plant residue com- 
partment (Figure 1), while dead wood goes into wood residue. Wood 
and plant residue decompose as a function of the lignin and N con- 
tent of the materia] and the abiotic decomposition pEirameter. The 
abiotic decomposition rate is a function of soil temperature and 
rainfall, and the equations presented in the Century model are used 
for this model (Figure 3). The abiotic decomposition rate increases 
with increasing temperature and increasing values of the monthly 
ratio of rainfall to potential evapotranspiration rate. The effect of 
plant lignin and N content on decomposition is represented by the 
functions used in the Century model. Plant and wood residue 
decomposition rates decrease with increasing lignin content and 
decreasing N content. As plant residue decomposes, most (60 to 
80%) of the C is lost to the atmosphere as respiration. The C lost as 
respiration is a function of soil texture and mineralogy and the plant 
lignin content of the vegetation. The effect of soil texture and lignin 
content on C0 2 respiration loss is represented using the equations 
presented in the Century model. C0 2 losses are higher for sandy 
soils and decrease as the lignin content of the vegetation decreases. 
Carbon not lost as C0 2 respiration flows into the two soil organic 
pools (slow and passive). The slow pool receives most of the carbon 
(98 to 99%). The fraction that goes to the passive pool is a function 
of the clay content of the soil, with higher values for clay soils. The 
slow and passive soil pools decompose at different rates (approxi- 
mate turnover times of 20 and 1000 years) which are pool-specific 
and vary as a function of the abiotic decomposition rate (Figure 3). 


‘Effect of ‘ Temperature (Tff 


308 


Molding the ‘Earth System 




Figure 3. (a) The effect of soil temperature (Jd) on monthly deposition rates, (b) 
The effect of moisture (M d ) on monthly deposition rates. 
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When the soil organic pools decompose, the carbon is returned to 
the atmosphere as C0 2 . 

Nutrient cycling is modeled by using the same flow diagram as in 
the carbon model. The nutrient content of the vegetation is an 
ecosystem-specific input variable, and the plant nutrients flow along 
with the carbon. When soil C pools decompose they release C to the 
atmosphere, and nutrients associated with the decomposing soil C 
are assumed to be mineralized and flow into available mineral pools 
(nitrate, ammonium, and phosphate: NO 3 -, NH 4 + and P0 4 s ). 

When plant residue decomposes it immobilizes nutrients from the 
mineral pools at a rate that is a function of the nutrient and lignin 
content of the vegetation. The immobilization rate is high er for vegeta- 
tion containing high lignin and low nutrients. The soil nutrient miner- 
alization model uses the equations developed in the Century soil 
organic matter model. The nutrients remaining after immobilization 
into decomposing plant residue are then available for uptake into new 
plant production. Nutrients will be lost from the system by leaching, 
by erosion, and through gaseous emissions. Nutrient input to the sys- 
tem comes from the atmosphere and N fixation by the plants. Simple 
equations for these Inputs and losses are included in the model. 

A representative description of soils for the present modeling 
effort and for interactions with atmospheric models will need to 
specify both hydrologic and fertility parameters. These will likely 
include depth, texture, presence or absence of an impedance layer 
(for water or root penetration), active carbon content, sum of 
exchangeable bases, and an additional description of topography. 

A combination of basic soil attributes such as parent material 
(from surface geological maps) and soil age, together with plant 
functional groups that grow on these soils (see below), may allow us 
to deduce both fertility and the hydrological cycle in an ecosystem 
without complex and detailed mapping of soil texture depth, infiltra- 
tion, interflow, and runoff. 


Biomes 

The ability to model different processes in the carbon and nutri- 
ent cycles of terrestrial ecosystems needs to be tested in different 
biomes. On a precipitation and temperature plot, eight major bio- 
mes can be distinguished, as listed in the previous section. Within 
each of these, on a regional basis (i.e., by general circulation model 
grid cell), it is necessary to define the proportion that has been con- 
verted to dryland or irrigated arable land. Feasible land use systems 
can be superimposed upon these biomes, and the model should be 
tested under these combinations. 
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Since the C model describes productivity in several key plant 
components, and also allows for the incorporation of nutrient fluxes, 
it can be used with only minor additional algorithms to predict other 
important ecosystem variables such as leaf area index, albedo, and 
to some extent moisture. 

A critical gap in this, as in most other models, is the inadequate 
representation of water use and storage within the ecosystem. 
Incoming precipitation controls primary production, but soil mois- 
ture storage and water uptake by roots from different soil layers at 
different times are not represented. The redistribution of precipita- 
tion into evaporation, soil moisture storage, transpiration, and 
drainage depends on such factors as relief, soil type, texture, rooting 
depth, and plant type. In sandy and loamy soils water will gradually 
penetrate the soil as successive layers become saturated. Water 
extraction (and loss to the atmosphere) subsequently depend on sur- 
face evaporation and root distribution. This Is relatively simple to 
represent in a model. In heavier soils, cracks permit the rapid infil- 
tration of incoming rainwater to depth without the complete satura- 
tion of upper soil layers. Subsequent water uptake and transpiration 
by plants from lower soil layers is an important factor for plant sur- 
vival in dry climates, but is difficult to represent with simple models. 
The complex redistribution patterns of water on the Individual 
ecosystem level will have to be simplified to provide water balances at 
a resolution appropriate at the general circulation model (GCM) grid 
scale. No satisfactory solution to this problem is yet available. 


Scenarios 

The utility of the carbon flow model for predicting C dynamics 
under global change depends on whether it can be applied to differ- 
ent scenarios of global change around the globe. The scenarios 
include different conditions and changes in atmospheric C0 2 levels; 
temperature; precipitation; atmospheric deposition of N, S, and oth- 
ers; land use, such as fertilization, irrigation, and cultivation; defor- 
estation/afforestation; extraction of selected ecosystem products; 
grazing; and fire. 

C0 2 levels are covered in the plant production portion oi the 
model structure. In Running’s model, photosynthesis depends on 
COo gradients that inherently depend on atmospheric C0 2 levels. In 
the Century model, the slope of the relationship between available 
water and productivity is modified by atmospheric C0 2 levels. 

Temperature and precipitation are controlling factors for both 
production and decomposition functions of the simplified Century 
models. In addition, the allocation of production to root and shoot 
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is controlled by moisture availability. Irrigation will have similar 
effects on production and decomposition through improved mois- 
ture availability. 

Atmospheric deposition may affect the availabilities of nutrients 
for plant production. Productivity is modeled depending on Cinutri- 
ent ratios, which allows for the effects of such atmospheric deposi- 
tion. Atmospheric deposition may also change acidity levels in the 
soil depending on the soil's buffering capacity. The potential produc- 
tivity depression will lower the slope of the curve relating productiv- 
ity to moisture availability of the simplified Century. 

Levels of inorganic nutrients in the soil as affected by atmospheric 
deposition or land management (including fertilization) may need to 
be modeled separately from the Cmutrient function. Phosphorus, for 
Instance, has several inorganic forms of differing solubility and avail- 
ability, which interact with organic P forms and plant uptake. Trans- 
formations between these forms need to be modeled separately. This 
has been done successfully in the complete version of Century. 

The effect of cultivation is to mix residues into the soil and speed 
up decomposition processes. This is successfully modeled by Cen- 
tury, and adequately represented in the present model. 

The processes of deforestation and afforestation Involve depletion 
and buildup of leaf, wood, and root components, which can be ade- 
quately modeled. 

Extraction of selected products such as firewood, fruit, or animals 
in grazed systems without wholesale modifrcation of the ecosystem 
can be represented by loss functions on the leaf or wood compo- 
nents. The partial return of materials through feces of grazing ani- 
mals or accelerated litter production in extractive use may require 
additional simple models that provide new parameters to the simpli- 
fied Century. 

Similarly, fire causes a redistribution of model components, 
export of C and N, and generation of inorganic forms of nutrients in 
ashes that will need to be modeled in a separate fire model, which 
should account for the differential effects of bum intensities, fuel 
availability, and vegetation state. 


Vegetation Change 

The eight biomes listed above may be exchanged for a more 
appropriate classification of vegetation functional types once these 
have been developed using plant functional types (PFTs) and soil 
properties. The PFTs will be defined in terms of characteristics that 
determine their interaction with the atmosphere and other ecosys- 
tem components. 
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Thus far the terrestrial carbon model discussed here Is applicable 
to a fixed vegetation cover within each of the eight blomes^ As climate 
Sid atmospheric C0 2 change, however, the functioning of the tares- 
trial ecosystem will change, and this will eventually lead to global 
changes in the state of the ecosystem. To take these changes into 
account, fine-scale models of detailed changes in species composition 
are inappropriate. What is needed are a broad-scale data base of g 
referenced characteristics of the earth's vegetation cover and a model 
that predicts changes that are significant in terms of their conse- 
quences for GCMs and atmospheric chemistry (i.e., carbon). 

We propose the development of two complementary models for the 
changing ecosystem: a state-and-transition model and an environ- 
mentally driven global vegetation model. 

State-and-Transition Model 

This model should be based on a present-day land cover map of 
the world, classified into the eight biome types, with local and 
regional differences represented by: 

. Separate descriptions of the vegetation functional types ^(wooch 
land, grassland, etc.), with each vegetation type being defined y 
its proportional composition of plant functional types (evergreen 
vs. deciduous, annual vs. perennial, etc.) 

• Differences in the amounts of the variables in the carbon model 
. The proportion of the biome that has been converted to cultivated 
land. 

The first requirement is that the model must apply to the entire 
land surface of the globe, and this puts a constraint on the scale^ 
grid cell or polygon size of about 100 x 100 km is tentatively pro- 
posed. For each grid (polygon) cell, the proportion of each biome 
type should be determined. It is necessary to determine the possible 
states in which each biome type can exist. The definition of these 
states should be guided by the interaction of the ecosystem with th 
atmosphere. Minor differences in species composition should not be 
considered as differences in state. Large differences in the propor- 
tions of different functional types of species (e.g.. from evergreen to 
deciduous trees, or from trees to grass), which will lead * 
able differences in the atmosphere/climate models, are signifleanL 
Transition from one state to another is driven by both endoge 
nous (community) processes and exogenous (disturbance) e ec 
Both need to be taken into account. Change due to endogenous 
processes will be determined mostly as a function of the changes 
feveTs of the C boxes and fluxes in the C model, but may also be 
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based on known successional processes in that type of ecosystem. 
In terms of environmental (precipitation and temperature) changes, 
present-day correlations between vegetation state and environment 
may be used to guide the degree of change In environment required 
to induce a “significant" change in state (with due recognition given 
to the various lag effects). In addition, there should be an attempt to 
derive relationships (within each blome) between the levels of C in 
each box of the model and the “state" of the vegetation. These rela- 
tionships could then be used to drive the vegetation change model. 

For ecosystems where change is predominantly a function of dis- 
turbance regime, the combinations and sequences of conditions for 
each transition need to be explicitly stated, based on what is cur- 
rently known about the ecology of that region. 

The two types of processes are then amalgamated into a set of 
rules which determines each transition. 

Environmentally Driven Global Vegetation Model 

We envisage this model as based on first principles, using funda- 
mental processes of vegetation development, and having no recourse 
to existing maps. It is based on a set of plant functional types with 
defined environmental responses, and the proportional composition 
of the vegetation is determined by the individual and interactive 
dynamics of these PFTs in response to given climatic scenarios. The 
proposed model is a simplified FORET-type model (Shugart, 1984) of 
plant-by-plant replacement, with globally defined plant functional 
types replacing individual species. 

Validation 

Given our present understanding of global vegetation dynamics, 
and the existing data (or lack thereof), it will be particularly valuable 
to develop both of these vegetation change models. Discrepancies 
between them will focus attention on high-priority areas for reduc- 
ing uncertainty. The vegetation model will be particularly valuable 
in setting bounds for the state-and-transition model and for check- 
ing on climate-driven changes to vegetation/environment combina- 
tions for which there are no current analogs. Other means of valida- 
tion will include the use of advanced very high resolution radiometer 
(AVHRR) data and existing data sets on vegetation, soil, and climate 
for various regions within each biome. 
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